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Abstract 



The E885 collaboration utilized the 1.8 GeV/c K~ beam line at the Alter- 
nating Gradient Synchrotron (AGS) to accumulate greater than 10 times the 
world's existing data sample of (K~,K + ) events on carbon. A total of about 
3 x 10 5 (K~,K + ) events were collected and analyzed. S hypernuclear states 
are expected to be produced through the reaction K~ + 12 C — > K + +^ Be. 
A signal could also result from direct production of j^Be + A without a dis- 
tinct E intermediate state. The measured missing-mass spectrum indicates 
the existence of a signal below the threshold for free E~ production that can- 
not be explained by background or effects of limited resolution. Although 
the resolution was not sufficient to resolve discrete hypernuclear states, the 
excess of events in the region of missing mass, kinematically inaccessible in 
free E~ production, can be compared with theoretical predictions for ^Be 
production. Reasonable agreement between the data and theory is achieved 
by assuming a H-nucleus potential well depth Vos of about 14 MeV within 
the Woods-Saxon prescription. 

PACS number(s): 21.80.+a, 21.30.Fe 
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I. INTRODUCTION 



With the advent of magnetic spectrometers, it became possible to study hypernuclei 
through the reconstruction of missing mass. Hypernuclei with one strange quark (S=-l 
sector) have long been studied using this method quite successfully. Similarly to production 
of S=-l hypernuclei in (tt + ,K + ) and (K~,7r~) reactions, the (K~,K + ) reaction can be used 
to create S=-2 hypernuclei. 

Dover and Gal found the S-nucleus potential well depth to be ~ 21 MeV - 24 MeV based 
on their analysis of emulsion data [ij. The Is and lp single-particle S~ states are expected to 
be bound for the values of the potential within that range. Although bound S-nucleus states 
would be broadened by the conversion EN — ► AA, calculations have indicated that their 
widths could be a few MeV or less M. Fukuda et al. recently reported evidence for bound 
5 hypernuclei produced in the (K~,K + ) reaction at px- = 1.6 GeV/c on a scintillating-fiber 
target (KEK E224) [^J. Their missing- mass resolution was not sufficient to resolve discrete 
hypernuclear states but their data seemed to favor a S-nucleus potential well depth of ~ 
16 MeV. They concluded that a value of 24 MeV for the potential was improbable. 

The 1.8 GeV/c K~ beam at the AGS was used to create doubly strange (S=-2) systems 
in E885 via the 12 C(K~,K + )X reaction at p K - = 1.8 GeV/c. A high K~ flux and a thick 
diamond target || combined to give a sensitivity to S=-2 states of 1.5 counts per nb/sr as 
compared to 0.05 counts per nb/sr for E224 ||. Although the experimental resolution in 
E885 was somewhat better than the E224 resolution (4.2 MeV/c 2 r.m.s. vs. 5.6 MeV/c 2 
r.m.s. for hydrogen kinematics), it too was insufficient to resolve discrete hypernuclear 
states. Nevertheless, an estimate of the cross section can be made; the excitation-energy 
spectrum can be compared to the theoretical predictions for different potential well depths 
folded with our experimental resolution. [] 

The organization of this paper is as follows. In Sec. [Tl] the experimental setup for E885 
is described and the event reconstruction is briefly discussed. In order to make cross sec- 
tion estimates and compare the data with theory, the missing-mass spectrum needs to be 
normalized. This normalization procedure is described in Sec. 111. Prior to a comparison 
with our data, the theoretical cross section calculation must be folded with the experimental 
resolution function. The determination of the resolution function is described in Sec. |TV]. 
The theoretical cross sections for a series of values of the S-nucleus potential well depth 
parameter, folded with the experimental resolution function, are compared with the data in 
Sec. [V|- A discussion of background is presented in Sec. [Vl|, followed by our conclusions in 

Sec. ivrq . 



1 Excitation energy E is defined as the missing mass minus the combined mass of the 11 B(gs) and 
the H~ (all times c 2 ). With this definition the threshold for free H~ production is at zero excitation 
energy and bound H hypernuclear states have negative excitation energy. 
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II. EXPERIMENTAL SETUP 



A side view of the E885 detector apparatus is shown in Fig. [IJ The K~ beam was delivered 
by the AGS D6 beam line 0] at a rate of 10 6 K~ /spill with a 1:1 K/-7T ratio. The spill duration 
and repetition rate were 1.6 s and 1000 spills/hour, respectively (44% duty factor). The 
purity and intensity of the kaon beam provided by the AGS D6 beam line greatly exceeded 
that available elsewhere. The beam was focused in the vertical direction by the last beam line 
quadrupole magnet onto a synthetic diamond target, 8 cm wide x 1cm high x 5 cm thick 
with a density of 3.3 g/cm 3 . Diamond was used as the target material in order to enhance the 
H~ stopping rate (a property not important for the direct production discussed in this paper); 
the resulting compactness of the target also allowed the application of tighter geometric cuts. 

The momentum of the incident particle was measured using information from a scintil- 
lator hodoscope, MP, located in the beam line at the exit of the first mass slit (not shown 
in Fig. [I|) and information from 3 drift chambers (ID1-3) located downstream of the last 
beam line magnet. The amplitude in the IC Cherenkov counter, along with the measured 
particle momentum and the beam line time-of-flight, were used to suppress triggers from 
pions. Additional details regarding the experimental apparatus can be found in Ref . || . 

The momentum of the outgoing K + was measured in the 1.4 T 48D48 magnetic spectrom- 
eter. The amplitude in the FC Cherenkov counter, along with the measured momentum and 
IT-BT time-of-flight, was used to select kaons. In addition to pions, protons were another 
source of background in the spectrometer. To reduce the hardware trigger rate, an aerogel 
Cherenkov counter, FCO, was used to suppress protons; its refraction index was selected to 
be sensitive to kaons but not protons from inelastic events. The measured mass spectrum for 
a sample of secondary particles is shown in Fig. [| This shows that the kaons were cleanly 
separated from the pions and protons that passed the hardware triggers after goodness-of-fit 
and vertex cuts were applied. The pion peak in Fig. 0a is due to leakage of the hardware 
FC pion veto. These events were virtually eliminated by an additional software cut on FC 
pulse height. The results are shown in Fig. 0b. The measured momenta of the incident and 
secondary particles allow the calculation of the missing mass. 



III. NORMALIZATION 

Our data were normalized to the known forward-angle cross section of free S~ production 
on hydrogen (contained in a CH 2 target). Due to the low statistics in the bound region, 
only angle-averaged differential cross sections could be obtained. These angle-averaged cross 

sections are denoted here as ( — — — ). We present results averaged over 0° < 6 K + < 14°, 

\ dildE I 

corresponding to the full acceptance of the spectrometer, and for 0° < 9 K + < 8°. 

The normalization effectively scaled the diamond target results to agree with the S~ 
production off carbon in a CH2 target, which was, in turn, normalized to the S production 
off hydrogen in the CH2 target. This procedure corrected for differences in the event recon- 
struction efficiency between the low-intensity CH 2 running and the high intensity diamond 
target running. These differences were brought about, in part, by tight multiplicity cuts on 
detector MP. 
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The scale for the excitation-energy spectrum was set using the equation: 



/ d 2 a \° _ fdN\ C 

\dndE/ ~ \lE) i - (1) 



diam 



where the normalization factor is given by 

A -- 



<* H Ag H2 NS H2 lda\ H 



The following definitions were used: 



/ d 2 a \° 



angle-averaged double-differential cross section for hypernuclear produc- 



\dVLdE I 
tion on carbon. 

( dN\ ^ 

— — - number of hyp ernuclear-pro duct ion events in an excitation energy inter- 
ns;, 



diam 

val dE from the diamond target. 

• a H - angle- and energy-averaged acceptance for free S~ production on hydrogen. 

• a c (E) - angle- averaged acceptance for hypernuclear production on carbon. The ar- 
gument E indicates the dependence on excitation energy. 

• ^ch 2 ~ target thickness for hydrogen in CH 2 target. 

• Xq H - target thickness for carbon in CH 2 target. 

• Nq H - number of events of free 2~ production on carbon in CH 2 target. 

• N CH - number of events of free S~ production on hydrogen in CH 2 target. 

• Nj iam - number of events of free H~ production on carbon in diamond target. 
I da \ H 

• ( -— ) - angle-averaged differential cross section of S~ production on hydrogen. 
\du J 

Since the spectrometer acceptance for hypernuclear production varies very little with 
energy over the small region of interest about zero excitation energy, instead of the energy- 
dependent acceptance a c (E), the acceptance was calculated for a fixed excitation energy of 
— 10 MeV. The angle-dependent acceptance used to determine a (E) was generated with a 
Monte Carlo simulation and is shown as the solid line in Fig. [3]. 

A summary of the parameters used for normalization in Eq. (j^) is given in Table | for 
two different acceptances; the first one is full spectrometer acceptance {6k+ < 14°) and the 
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second is for limited acceptance < 8°) Q. The following values for the normalization 

factor A, giving the angle-averaged differential cross section per event were obtained: 



This factor relates the counts per MeV to the angle-averaged differential cross section per 
MeV. The statistical errors given for the normalizing factor A were due mostly to the limited 
amount of calibration data, namely the number of free H~ production events on hydrogen 
and carbon in the CH 2 target. The systematic uncertainty, excluding the uncertainty of 
the elementary cross section of 5~ production on hydrogen, is expected to be below the 
statistical error. The uncertainty of the elementary cross section for production on 
hydrogen does not affect the comparison between data and theory because the theoretical 
DWIA calculations are normalized to the same elementary cross section. Thus the data-to- 
theory ratio is insensitive to the assumed elementary cross-section. 



Data sets were taken with two different CH 2 targets and were used to study the missing- 
mass resolution. The first set used an 8 cm long target and the second set, which was 
acquired over a shorter time period, used a 13 cm long target. The first data set was used to 
study the spectrometer's intrinsic resolution and, in particular, to evaluate the tails of the 
resolution function. The second set was used to derive the width of the resolution function. 
The contribution of energy-loss effects to the resolution for the 13 cm long CH 2 target is 
expected to be similar to the contribution from the 5 cm long diamond target. 

Due to the two-body nature of H~ production on hydrogen, the missing mass spectrum 
is a peak, whose width equals the experimental resolution, on top of a background from H~ 
production on carbon. In addition to the peak width, we are also interested in the tail of the 
resolution. The study of the resolution tails is important; the small hypernuclear-production 
signal is not well separated in missing mass from the large signal of free S~ production. The 
background under the hydrogen peak, consisting of H~ production on carbon, is subtracted 
using the missing-mass spectrum of the carbon target. The resulting background-subtracted 
spectrum is shown in Fig. |j. The width of the peak from a Gaussian fit is 3.8 MeV/c 2 r.m.s.. 
It can be seen from the plot that the resolution is reasonably close to being Gaussian. The 
missing mass resolution for the 13 cm CH 2 target was measured to be 4.2 MeV/c 2 r.m.s. for 
hydrogen kinematics. 

The mass resolution is a function of the target mass. The resolution specific to 2~ 
production on hydrogen (which we measure directly) was used to infer the missing-mass 
resolution for carbon kinematics using simple kinematical considerations. A factor of 1.45 



2 Note that in both cases we use the value of 35 /ib/sr for the elementary production cross section 
which is compatible with the measured value in Ref. R. This factor should be divided out when 
comparing these results to a calculation of the effective proton number. 



.4 



( 



0.64 ± 0.05 nb/sr for 6 K + < 14° 
2.14 ± 0.20 nb/sr for 6 K+ < 8°. 



IV. MISSING-MASS RESOLUTION 
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relates the missing-mass resolutions in carbon and hydrogen kinematics. This value results 
in a missing-mass resolution for carbon kinematics of 4.2 MeV/c 2 x 1.45 = 6.1 MeV/c 2 . The 
peak of H~ production on hydrogen was also used for the energy-scale calibration. The 
position of the peak was adjusted to be at the S~ mass value and the accuracy of the 
energy-scale calibration for the diamond target was estimated to be better than 0.5 MeV. 



V. COMPARISON OF THE DATA AND THEORY 

The theoretical (K~,K + ) double differential cross section for ^Be production, , 
was calculated, for both bound and unbound S~s, in the Kapur-Peirls framework which was 
developed originally for an estimate of the (ir + ,K + ) strength function ||. The DWIA cal- 
culation was performed for a K~ momentum of 1.8 GeV/c by assuming the K-N elementary 
cross sections: a(K~ iV)=2.90 fm 2 and a(K + N)=1.99 fm 2 . Further details of the calculation 
can be found in the recent references and [Q in which the (K~,K + ) cross sections for the 
12 C and 16 targets have been estimated for p K - = 1.6 GeV/c. 

The calculation was performed for 0° < 9 K + < 16° and a series of Woods-Saxon H~ 
well depth parameter values: Voh = 12, 14, 16, 18, and 20 MeV. The radius and skin 
depth of the 5 potential were fixed at R = 1.L4 1 / 3 fm and a==0.65 fm. The proton wave 
function in the 12 C target was generated using a Woods-Saxon potential with V at = 50 MeV, 
R = 1.1A 1 / 3 fm and ajy=0.65 fm. The elementary H~ production cross section was set to 
35yub/sr to be compatible with the normalization of our experimental data. The kinematic 
factor a, which accounts for the transformation between the two-body and A-body frames, 
was set to 0.73. 

As a typical example of the angular dependence of the differential cross section, the 
dashed line in Fig. |3] shows the case of the ground state for Voh=14 MeV. In order to make 
a comparison between the theory and the data, we calculate the angle-averaged differential 

cross sections, ( — (E) ) and fold the results with the experimental resolution. However, 

\dndE K 7 F 

we first present theoretical 12 C(K~,K + )4 2 Be spectra which have not been folded by the 

experimental energy resolution but have been angle-averaged over 0° < 9k+ < 14°. As 

shown in Fig. [5], the result for Vo==20 MeV (dashed line) has two bound-state peaks, 

corresponding to the S~ s- and p-orbitals. The widths of these peaks are determined using 

a one-boson-exchange model to estimate the rate for the SiV — > AA conversion. For the case 

of Voh =14 MeV, the p state is not bound but it is calculated as a resonance state in the 

continuum; therefore a sudden rise is seen just above the threshold in Fig. [5]. 

Figure |6| shows experimental excitation energy histograms for 12 C(K~,K + )X for two 
different limits on the scattering angle of the outgoing K + , K + < 14° and 6k+ < 8°. The 
data clearly show an enhancement around zero excitation energy when compared to a Monte 
Carlo simulation based on quasi-free H production which has been normalized to the total 
number of 12 C(K~,K + )X events (curve QF). 

In the same figure, the ^Be production theoretical curves for several H potential well 
depths, folded with the 6.1 MeV r.m.s. experimental resolution, are shown for comparison 
with the data. The expected location of the ground state of A \Be (assuming a total binding 
energy of the A's, Baa, of 25 MeV) and the thresholds for A x Be + A and n B + H~ production 
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are indicated. The normalization calculation for the case 6k+ < 8° is less sensitive to the 
model of angular dependence because the spectrometer acceptance is fairly flat over this 
region but drops rapidly for 9%+ > 8° as shown in Fig. |3|; we present the results for both 
the entire acceptance and for 9 K + < 8° in Fig. ||. 

Visual inspection shows that the theoretical curve for the value of the S-nucleus potential 
well depth Vos = 14 MeV agrees with the data reasonably well in the region of excitation 
energy -20 MeV<E<0 MeV and much better than the curve for Vo~ = 20 MeV. If any 
of the observed signal results from direct two-A production without an intermediate i 2 Be 
state, the discrepancy between the Vo= = 20 MeV results and the remaining experimental 
signal becomes even larger. The comparison of data to theory in the bound region can be 
quantified by considering the results for the angle-averaged cross section integrated over the 
-20 MeV<E<0 MeV excitation region. The results are presented in Table [TI|. It can be 
seen that the DWIA calculation predicts a yield more than a factor of two larger than the 
measured cross section in this region when a S-nucleus potential well depth of 20 MeV is 
used. 

VI. BACKGROUND 

Events due to S production off the hydrogen contained in the target and target-area 
scintillators produce a small peak at ^+90 MeV in excitation energy; this peak represents 
only a few per cent of the total events and thus this contamination does not affect the analysis 
near zero excitation energy. Background in the region of negative values of excitation energy 
can be divided into two groups: 

• Flat background caused by particle misidentification, severe tracking errors, etc. Such 
background might have some structure over a wider region of missing mass, but on 
a small, ~ 100 MeV scale, it should be flat. This was confirmed when observing 
the background in the bound region obtained by using a loose set of cuts (hence 
considerable background contamination). 

• Leakage from E>0 due to finite missing-mass resolution. 

Background of the first kind can be easily estimated by observing the region further below 
zero where leakage becomes negligible. There are 6 events in the interval of excitation 
energy [-80 MeV,-40 MeV] for 9 K + < 14° (see Fig. Since the flatness of the background 
was confirmed by loosening the cuts, an extrapolation can be made that the number of 
background events of the first kind in the excitation-energy region [-20 MeV,0 MeV] is 
about 3. It can be concluded then that the background of the first kind contributes only a 
small fraction of the total signal in the S hypernuclear production region, and it is ignored 
when comparing the data to theory. 

Background of the second kind is taken into account when the theoretical cross section 
is folded with the experimental resolution. The number of events in the E<0 region that 
leak from the E>0 region depends on the shape and magnitude of the excitation-energy 
spectrum for E>0 and on the resolution function shape. We found that the signal observed 
in the bound region cannot be explained by leakage from the E>0 region. This is already 
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evidenced by the good agreement between the data and theory for Vo= = 14 MeV in the 
E<0 region, since non-zero cross section of hypernuclear production is essential to explain 
the events in the bound region. We, however, made two additional checks to demonstrate it. 
First, a Monte Carlo excitation-energy spectrum, where only free S~ production is simulated 
and no events in the E<0 region are produced, was folded with the measured experimental 
resolution. It was found that only « 10 events are expected to have E<0 for 9k+ < 14° (the 
observed number is 67). Second, a different resolution function was tried with very large 
tails added to the basic Gaussian resolution function. Even with those large tails, which 
are clearly ruled out by the resolution measurement, the signal in the bound region due to 
leakage was still a few times below the observed level. 



VII. CONCLUSIONS 

The measurement of the missing mass in the (K~,K + ) reaction on carbon in E885 allowed 
the determination of the excitation-energy spectrum of directly produced S=-2 systems. The 
signal in the missing-mass region below the threshold for free H~ production was examined 
in order to extract information on the existence of S hypernuclei and the strength of the 
S-nucleus potential. Our analysis has shown that the events in the bound region of the 
missing-mass spectrum could not be explained by non-S=-2 background or effects of limited 
resolution, and these events are consistent with hypernuclear production. Possible sources of 
the signal in the bound region are the production of one or two As in the continuum or double- 
A hyperfragments. The threshold for j^Be + A is estimated to be at E « -27 MeV. This 
channel may dominate the other two-A channels since the two-step (K~,7r),(7r,K + ) process 
involves one low-momentum A and one high-momentum A. Its cross section is expected to be 
several orders of magnitude larger than the production cross section for double-A hypernuclei 
10|,11|. No quantitative estimate of the contribution of these processes in this kinematic 



regime has been made and final interpretation of the excitation spectrum in terms of one- 
step S-hypernuclear production verses two-A production without an intermediate S-state 
requires additional theoretical work. 

A comparison of the data and the theoretical calculations for the cross sections of H hyper- 
nuclear production, folded with the experimental resolution, shows a reasonable agreement of 
the data and theory for Vos = 14 MeV in the region of excitation energy -20 MeV<E<0 MeV 
and a significant disagreement between the data and theory for Voh = 20 MeV. These re- 
sults are consistent with the conclusions of Ref. H]. If a significant portion of the events 
in this region correspond to two-A production without an intermediate 5 hypernucleus, the 
discrepancy between the calculation using this well depth and our results would be even 
larger. We conclude that the results are consistent with the theoretical predictions when a 
potential depth of 14 MeV or less is assumed. 

The accuracy of DWIA, which was used for the cross section calculation, was nominally 
assumed to be 30% in Ref 0. Although the DWIA framework has been successful in 
describing the production of single-A hypernuclei using the (7r + ,K + ) and (K~,7r~) reactions, 
it is possible that that the DWIA accuracy for the (K~,K + ) reaction is worse due to the 
higher momentum transfer. Therefore a value for the H-nucleus potential well depth as high 
as 20 MeV, though less favored, cannot be completely ruled out. Additional theoretical 
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calculations of the complete 12 C(K~, K + )X excitation spectrum are needed to determine the 
parameter space that is compatible with our results. 
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TABLES 



maximum 6^+ 
(degrees) 


a 11 


a c (-10 MeV) 


A CH 2 1 A CH 2 




(da/dn) H 
Oub/sr) 


^Eh 2 


14 


8.1% 


14.8% 


2 


1/125.6 


35 


480 


8 


14.4% 


15.3% 


2 


1/127.9 


35 


240 



TABLE I. A summary of parameters used for normalization. The value used for the cross 
section of S~ production on hydrogen is compatible with the value measured by T. Iijima et al. 0. 



This Work KEK E224 





e K + <8° 


e K+ <i4° 




Counts 


42 


67 


3 


Sensitivity (counts per nb/sr) 


0.47 


1.6 


0.05 


(da/dn) Experiment (nb/sr) 


89±14 


42±5 


60±45 


(da/dn) V 0S =12 MeV (nb/sr) 


67 


26 




(dff/dtl) V 0S =14 MeV (nb/sr) 


92 


37 




(da/dn) V 0S =16 MeV (nb/sr) 


115 


46 


75 


(da/dn) V 0S =18 MeV (nb/sr) 


166 


68 




(da/dn) V ==20 MeV (nb/sr) 


226 


93 


170 



TABLE II. The experimental results integrated over the excitation energy -20 MeV<E<0 MeV 
are shown. The theoretical values for the averaged differential cross sections, after convoluting with 
the experimental resolutions and integrating over the same energy range, are shown as a function 
of the H-nucleus potential well depth parameter, Vos- The corresponding values for the KEK E224 
experiment ||, p^-=1.66 GeV/c, are shown where available. 
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FIGURES 




Target 48D48 



-1 01 2345678 
I I I I I I I I I I 

Z(m) 

FIG. 1. Detector configuration for E885. The drift chambers ID1-3 determine the incident K~ 
trajectory and, combined with the beam line hodoscope MP (not shown here), the K~ momentum. 
The drift chambers FDO-3 and BD1-2 determine the K + trajectory through the 1.4 T 48D48 dipole. 
Scintillators IT and BT determine the K + time-of- flight. Hodoscopes FP and BP determine the 
spectrometer acceptance. Aerogel Cherenkovs IC, FC, and BC reject pions, while FCO rejects 
protons. Above and below the target are scintillating fiber detectors. On the left and right of 
the target are neutron detectors. The fiber detectors and neutron detectors were not used in this 
analysis. 
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0.25 0.50 0.75 1.00 

Secondary Particle Mass (GeV/c 2 ) 

FIG. 2. The secondary particle mass as calculated from the measured momentum, path length, 
and time-of-flight through the K + spectrometer for a subset of the data, (a) Events that pass 
goodness-of-fit, vertex, and multiplicity cuts, (b) An additional software veto has been applied 
using the FC pulse height. Events in (b) that fall within the secondary-mass cut were used in the 
analysis. 



(degrees) 
12 8 




0.00 



0.97 0.98 0.99 1.00 

COS 

FIG. 3. The angular acceptance, a(9), generated by a Monte Carlo simulation is shown as the 
solid line. The dashed line shows the results of the DWIA calculation of the angular dependence 
of the cross section of the ground state using a potential depth of Vos=14 MeV. 
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1.3 1.31 1.32 1.33 1.34 1.35 1.36 1.37 



Missing Mass (GeV / c 2 ) 

FIG. 4. Missing mass for H~ production on hydrogen in the 8 cm long CH2 target. Background 
from reactions on carbon has been subtracted using a data sample, properly normalized, taken with 
diamond target. The solid line shows the best Gaussian fit. 
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FIG. 5. Results of DWIA calculations, before folding by the experimental energy resolution, for 
the 12 C(K~, K + )i 2 Be reaction for Vo~=14 MeV and 20 MeV. The cross section has been averaged 
over the kaon angular range < 9 K + < 14°. 
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FIG. 6. Excitation-energy spectra from E885 for 12 C(K~,K + )X for 9^+ <14° (top figure) and 
K + <8° (bottom figure) along with ^Be production theoretical curves for Vos equal to 20, 18, 16, 
14, and 12 MeV. The results of a quasi- free S production calculation are also shown (curve QF). 
The expected location of the ground state of ^A^e anc ^ the thresholds for j^Be + A and 11 B + 
production are indicated with arrows. 
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